Jamming Transition and Inherent Structures of Hard Spheres and Discs 
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Recent studies show that volume fractions at the jamming transition of frictionless hard 
spheres and discs are not uniquely determined but exist over a continuous range. Motivated by this 
observation, we numerically investigate dependence of tpj on the initial configurations of the parent 
fluids equilibrated at a fraction ip cq , before compressing to generate a jammed packing. We find 
that <pj remains constant when (p oq is small but sharply increases when ip C q exceeds the dynamic 
transition point which the mode-coupling theory predicts. We carefully analyze configurational 
properties of both jammed packings and parent fluids and find that, while all jammed packings 
remain isostatic, the increase of ifij is accompanied with subtle but distinct changes of (i) local orders, 
(ii) a static length scale, and (iii) an exponent of the finite size scaling. These results quantitatively 
support the scenario of the random first order transition theory of the glass transition. 

PACS numbers: 64.70. P-,61.43.Fs, 45.70.Cc 



Despite of their apparent similarities, a unifying the- 
ory of the glass transition of supercooled fluids and the 
jamming transition of athcrmal particles such as granu- 
lar materials is still missing. Both are characterized by a 
transition from flowing state to randomly jammed state 
at a finite density or temperature. For frictionless par- 
ticle systems (which we shall consider in this Letter), it 
has long been argued that the jamming transition is in- 
terpreted as the zero-temperature limit of the glass tran- 
sition Numerical studies, however, show that the 
nature of two transitions is more complicated. For ex- 
ample, the jamming transition has been believed to take 
place sharply at an unique volume fraction in the ther- 
modynamic limit, the so-called "points J"; ipj ~ 64% 
for three dimension (3D) and 84% for two dimension 
(2D) Q. But recently it has been demonstrated that (pj 
is not unique but exists over a continuous range of volume 
fractions whose values vary depending on the protocols 
used to generate the jammed states [3-5]. The jammed 
configurations at different if], however, are found to re- 
main isostatic, lack a partial crystalline order, and there- 
fore are not mixtures of ordered and "maximally random 
jammed" (MRJ) states [111]. 

On the other hand, our understanding of the glass tran- 
sition is no better than that of jamming transition. Even 
a mean field picture of the glass transition has not been 
established. A promising candidate is the so-called ran- 
dom first order transition theory (RFOT), originalh 
spired by the mean-field theory of the spin-glasses 
Crudely speaking, the RFOT integrates the energy land- 
scape picture, the concept of the ideal glas s transition, 
and the mode-coupling theory (M CT) [§- 11 1 . Despite its 
theoretical coherence, the RFOT-MCT scenario still re- 
mains controversial, partly due to the lack of impeccable 
numerical and experimental evidences. 

The goal of this Letter is to provide numerical evi- 
dence that the protocol dependence of if] is a natural 
consequence of RFOT and MCT and thus the scenario 



can unify the glass and jamming transitions of friction- 
less particles. The idea that the free-energy landscape of 
glasses is intimately related to the jamming transition is 
not new 0, 12 - 14 1 . But to the best of our knowledge, 
quantitative characterizations of the transition points, 
particle configurations, and the associated length scales 
have not been done so far. First, we shall briefly re- 
capitulate the essence of RFOT-MCT scenario Q and 
how it integrates the jamming and glass transitions fl2l ]. 
In the mean field limit, RFOT predicts that, as a fluid 
is cooled down, it first undergoes the dynamical transi- 
tion at T mct followed by the thermodynamic transition at 
Tk(< T mct ). Below Tmct, the multidimensional (free) en- 
ergy surface becomes suddenly rugged. Concomitantly, 
properties of the energy surface abruptly change. The 
value of local minima of the energy surface or the in- 
herent structures (IS), e IS , which is almost constant at 
high temperatures start decreasing at T mct . At the same 
time, the saddles of the energy surface vanishes and all 
stationary points become stable. Dynamics near T mct is 
described by MCT, which predicts that dynamical quan- 
tities such as the relaxation time r Q diverge with a power 
law \T — T mct |~ 7 where 7 is a parameter also calculated 
by MCT [ll|. In finite dimensions, however, the dynamic 
transition is smeared out by activation hopping between 
local minima separated by finite barriers and becomes 
merely a crossover. Important observation is that the ge- 
ometrical properties of the energy landscape are not con- 
trolled by a single temperature T mct any more. Simula- 
tions have revealed that saddles survive well below T^£' , 
the dynamic transition point theoretically obtained from 
MCT, until they vanish at a lower temperature T t h, a 
so-called threshold temperature lB- 18|. The simulated 
relaxation time is still well fitted by MCT's power law, 
t q ~ \T— T^t\~' 1 but T^£t use d for fitting was found to 
be considerably lower than [ll| ■ Surprisingly, T, 
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are very close to T t h (TB - 17 1 . On the other hand, the on- 
set temperature T D at which e IS starts decreasing is close 
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to T, 
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[IE Oil. Downward departure of T^t (~ T th) 
from T^. e '(Rj T ), all of which should be identical in the 
mean-field limit, is due to the non-mean-field effect and 
can be explained using kinetic arguments [HED]. 
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The above argument also applies to hard sphere flu- 
ids. Temperature and the energy (T, e IS ), relevant vari- 
ables for continuous potential fluids, should be replaced 
by the (inverse) pressure P _1 and the volume for hard- 
core potential systems 22j. Instead of the volume, we 
shall adopt the density, or volume fraction ip. The inher- 
ent structures tp IS are obtained by compressing a parent 
fluid equilibrated at a finite P by letting P — > oo (with an 
extra minimization using a conjugate gradient method), 
just as T is quenched to zero to obtain e IS for the continu- 
ous potential fluids. This minimization process is nothing 
but generating the jammed packing for frictionlcss hard 
spheres and thus tp is should be equivalent with </?,j. Em- 
ploying the RFOT scenario discussed above, we predict 
that tp ls or if] is unchanged as long as P of the parent 
fluid equilibrated at tp cq is low but starts increasing as P 
(or ip eq ) exceeds P mc t (or ip mc t). In other words, ipj is not 
an unique value but is a function of P or ip cq and can ex- 



ists over a continuous range [l2-14|. The largest </>j max '' 
would correspond to the inherent structures of the fluid 
at the thermodynamic transition point, i.e., tp cq = tpx- 
For finite dimensional systems, like the continuous po- 



tential systems, <p 
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obtained from MCT theoretically 

(fit) 



should be lower than t^ct obtained by fitting the sim- 
ulation data for the relaxation times. Furthermore, tpj 
should increase if we prepare a dense parent fluid such 
that ip cq > V'mc? (rather than y>mct)- Another important 
prediction of RFOT is that at the dynamic transition 
point, the system enters to the coexisting region of nu- 
merous metastable phases, or mosaics. Thus, the static 
length scale associated with the mosaics, if any, should 
appear at ip^.^ ■ There have been several attempts to di- 
rectly measure the static length in supercooled fluids but 
most studies have focused on the configurations of parent 
fluids and at far lower temperatures or higher densities 
than the dynamic transition points Q. 

In order to verify these predictions, we prepare ther- 
mally equilibrated hard spheres (3D) and discs (2D) 
at various initial densities tp cq and study their inher- 
ent structures <pj. Both systems studied here are 50:50 
binary mixtures with a size ratio of 1.4 with periodic 
boundary conditions 0, 0] ■ The system is equilibrated at 
<p eq using Monte-Carlo simulation and then compressed 
rapidly to generate jammed states. Following the proce- 
dure employed in Refs. 0, (23[ , we switch the hard-core 
potential with the soft harmonic potential just before 
compression, allowing particles to overlap. The system 
is then relaxed to the zero-energy state using the con- 
jugate gradient method. This compression and energy- 
minimization cycle is iterated till the volume fraction is 
maximized without particle overlap. Note that the algo- 
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FIG. 1: (color online), <pj as a function of ip cq for binary 
mixtures of (a) 3D and (b) 2D. Arrows indicates the positions 

of V'mct, - Broken horizontal lines are ypj reported by O'Hern 
et al. 0]. 



rism to generate the jammed state is essential. For ex- 
ample, the Lubachevsky-Stillinger algorithm is inappro- 
priate as the system keep equilibrating during slow com- 
pression and find lower local minima of the landscape or 
higher tpj [24| . The system sizes are varied from N = 64 
to 2048. tpj in the large A^-limit is evaluated using the 
finite size scaling \ipj(N) — c/5j| ~ N~ x / vd . In Figure 1, 
the dependence of tpj on ip cq is shown. The exponents 
v = 0.72 (3D) and 0.74 (2D) obtained for the smallest 
ip cq are used for the rest of data. Actually, we found 



that v varies noticeably depending on ip cq as we shall 
discuss below. However, the different i/'s do not affect 
appreciably the results of FigQ] other than more scatter- 
ing of data points and larger error bars. At small ip eq , 
the jamming transition points are identical with those 
already reported in the literatures, v?j w 0.648 (3D) and 
ipj « 0.842 (2D) 0]. However, tpj abruptly starts increas- 
ing at large tp eq . The onset fractions are found to be very 

close to V'mct'* independently evaluated by solving the 
MCT equations for binary mixtures (indicated by arrows 
in FigHJ, ip { ^ ] « 0.516 (3D) and ^ t e) « 0.685 (2D). 
These values have been evaluated using the static struc- 
ture factor matrices S(k) obtained by the Percus-Yevick 
integral equation. Even if simulated S(k) are used, the 
values do not change more than 2%. The onset points 
are clearly much lower than y^ct ~ 0.59 (3D) and 0.79 
(2D) obtained from fitting the simulation data [2|| [2tl . 
Values of Fig.l are in accordance with those in Ref. [3j 
(3D) and consistent with Ref. [f| (2D). We also found 
that the nature of the jamming transition with larger 
<y9j's remains isostatic, i.e., the contact number at pj 
is given by z = 2d, where al is the spatial dimension, 
whereas the number of rattlers slightly increases as if] 
increases Q . We also measured the time sequence of the 
inherent structures for several <^ eq and observed that the 
patterns of sequences qualitatively change from whitc- 
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FIG. 2: Compositional order parameters (/ll, fss, /sl) for 
3d (a)~(c) and 2d (d)~(f). Vertical broken-dots lines repre- 
sent 95mct an d horizontal broken lines are guides for eyes. 



(not shown) , implying that the nature of the landscape 
is altered [27[. These results support quantitatively that 
the jamming and glass transitions can be discussed un- 
der the common rubric of the RFOT-MCT scenario and 
also that <p^£ is not a fictitious value of an approximate 
theory but bears the essential geometrical meaning. 

In order to clarify the nature of the denser jammed 
packings obtained from the parent fluid at (p cq > 
we focus on properties of their configurations. We calcu- 
late the compositional and orientational orders. Figured] 
shows the dependence on (p eq of the compositional order 
parameters of the jammed packings (/ll, fss, fsh), the 
number fraction of the contact pairs of the large (L) and 
small (S) particles For ideally random configurations, 
/ll + fss ~ fsh ~ 0.5 holds. Indeed, this is the case for 
Minute but sharp changes are observed at 
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•£n\ct- Especially, /ll {fss) for 3D increases (decreases) 
by about 5%, whereas variation of /ls is far smaller. For 
2D, qualitatively similar changes are observed, consistent 
with Ref. |||. We next analyze the bond-orientational 
order (BOO) parameters Q\ and Qq (3D) and (2D) 
defined in Rcfs. (HI |28|. The BOO paramters evaluated 
for the large particles are shown in Figure [3] (a) ^(c). The 
results for the small particles show qualitatively similar 
behaviors although the variations are less pronounced. 
All results demonstrate that the BOO parameters are 
constant at <^ cq < but changes abruptly at ^ct''- 

One may want to aruge that synchronized change of tpj 
and compositional/orientational orders is due to the on- 
set of a partial crystalization or demixing and that the 
system trances a line connecting smoothly the MRJ pack- 
ing at the smallest </?j and the maximally ordered and 




FIG. 3: Bond-orientational order parameters for large par- 
ticles of jammed configurations, Qe (a) and Q4 (b) in 3D 
and for 2D (c). The inset shows ^6 for the parent fluid 
(Eq; empty squares) compared with the data for the IS (filled 
circles), (d) The correlation length ^6 obtained from the cor- 
relation function ge(r) of the jammed configurations (IS, filled 
circles) and of the parent fluid (Eq, empty squares) for 2D. 

Arrows represents f^f^ ant ^ horizontal lines are guides for 
eyes. 
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However, as demonstrated by 
the variations in both and 
the BOO paramters are far smaller than that expected 
from the coexisitence of partial order, suggesting that 
the system is riding on a different branch. Of course, 
we did not observe any hint of the local order from the 
eye inspection of the jammed configurations. We em- 
phasize that the sharp changes at ip^l are only ob- 
served in the inherent structures. The inset of Figj3] (c) 
shows that ^6 of the parent fluid continously incereases 
with ip cq with no hint to change around (p^^ . Simi- 
lar results were obtained for 3D. According to RFOT, 
the sudden decrease of the free-enegy at v^ct^ heralds 
the appearance of numeous metastable states and the 
system "phase-separate" into many amorphous states or 
mosaics. Thus, it is expected that these mosaics and their 
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FIG. 4: i/Peq-dependence of v of the finite-size scaling law 
\(pj(N)-ipj\ ~ Af" 1 ^ for 3D (a) and 2D (b). The horizontal 
broken lines are the values reported in Ref. 0]. 
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associated length scale should appear at (p 
tect such length, we calculate the static correlation func- 
tion of the fluctuations of the BOO parameters such as 
96(f) = (<5^ , 6( r ) ( 5 v f , 6(0)} and extract out the length scale 
^6 by fitting the results with the Orstein-Zernike func- 
tion. The reuslt for 2D is shown in Figure [3] (d). Similar 
result was obtained for 3D. ^6 which is almost unity at 
low (p e q starts incrasing at y^t^. Also shown is £g ob- 
tained for the parent fluid, which monotonically increases 
with tp cq . Although it is premature to conclude, these 
results suggest a posibility that the increase of ^6 of the 
jammed packing is a direct reflection of the emergence of 
the mosaics. 

Finally, we argue that V?mct ^ may also mark the point 
beyond which the finite-size scaling law is qualitatively 
altered due to the emergence of mosaics. It has been 
argued that in the crossover region at which the MCT's 
critical dynamics and activation processes coexist, the fi- 
nite size effect is highly nontrivial [2{|. For the short- 
range interaction systems, these two mcchanims may 
compete and the simple power-law scaling may be vio- 
lated. Figure |4] shows (^oq-dependence of the finite size 
scaling exponent v, obtained by naively using the scal- 
ing law. v is constant at ip eq < ^mc? and close to the 
value reported in Ref. Q, but v starts fluctuating and 
becomes errant at ip cq w <p^£ ■ We presume that this is 
another, though indirect, evidence supporing the RFOT- 
MCT scenario. 

In summary, we have accumulated and displayed 
qunatitative evidences that the RFOT-MCT scenario in- 
tegrates the jamming and glass transitinons in a com- 
mon language and successfully explains the continuous 
increase of </?j reported previously. We demonstrated 
for the first time that the dynamical transition point 

Vmc? theoretically evaluated, and not v?mct obtained by 
the fitting, unambiguously marks the onset of qualita- 
tive changes of the free energy landscape of hard spheres. 
Note that the results shown here are consitent with re- 
sults reported for short-ranged potential systems 11911 but 
not with those of the fully-connected models [30|. In 
Ref. [30l ] , the onset volume fraction at wihch tpj start in- 
creasing is considerably lower than <p^£ obtained from 



simulated relaxation time in the mean-field limit. This 
contradictory result might be attributed to the fully- 
connctcd intcarctions of the model. Indeed, it is known 
that the onset temeprare of the inherent structures for 
a fully-connected spin-glass model at finite sizes is much 
higher than the mean-field values and the convergence to 



the mean-field limit is extremely slow [31 



All results in this Letter eloquently support the RFOT- 
MCT sceniro but many nagging questions are left for us. 
For example, why does the MCT work quantitatively so 
well in finite dimensions? It is especially puzzling be- 
cause recent studies show that the traditional MCT is 
not perfectly consistent with the mean-field scenario at 
large spatial dimensions [HJ . Also we are left unanwered 
about the relation of the static length which we observed 
with other lengths observed via static or dynamic mea- 
surements in the past [H [H, 0|. And last interesting 
question may be whether the configurational properties, 
especially isostatic nature, of jammed packings is affected 
when (p cq exceeds (p^ll a t which all saddles of the free 
energy surface near the IS vanishes. These are a few of 
many problems which are left for the future works. 
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